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ABSTRACT 

V621Persei is a detached eclipsing binary in the open cluster xPersei which is com- 
posed of an early B-type giant star and a main sequence secondary component. From 
high-resolution spectroscopic observations and radial velocities from the literature, we 
determine the orbital period to be 25.5 days and the primary velocity semiamplitude 
to be K = 64.5 ± 0.4 km s -1 . No trace of the secondary star has been found in the 
spectrum. We solve the discovery light curves of this totally-eclipsing binary and find 
that the surface gravity of the secondary star is logge = 4.244 ± 0.054. We compare 
the absolute masses and radii of the two stars in the mass-radius diagram, for differ- 
ent possible values of the primary surface gravity, to the predictions of stellar models. 
We find that log<?A ~ 3.55, in agreement with values found from fitting Balmer lines 
with synthetic profiles. The expected masses of the two stars are 12 M and 6M and 
the expected radii are 10 Hq and 3 Pi©. The primary component is near the blue loop 
stage in its evolution. 

Key words: stars: binaries: eclipsing — open clusters — stars: fundamental param- 
eters — stars: binaries: spectroscopic — stars: distances — stars: early-type 



1 INTRODUCTION 



Detached eclipsing binaries are fundamental sources of data 
on the properties of stars (Andersen 1991). The masses and 
radii of two stars in a detached eclipsing binary (dEB) can 
be empirically determined to accuracies of better than 1% 
(e.g., Torres et al. 2000). Accurate effective temperatures 
and absolute magnitudes can be found using empirical pho- 
tometric calibrations, and the distances to dEBs can likewise 
be found to accuracies of 5% without any use of theoretical 
calculations (e.g., Clausen 2004). 

As the component stars of a dEB have the same age, 
chemical composition and distance, their physical properties 
provide an important test of the predictions of theoretical 
stellar evolutionary models. This test is aided if the stars are 
quite dissimilar in mass and radius, and if they have evolved 
away from the main sequence. Unfortunately, as the age and 
interior chemical abundances of stars are not, in general, 
observable, stellar models of any chemical composition and 
age can be compared to the observed data in order to find 
a good fit. 

Eclipsing binaries in open clusters provide a way of find- 
ing the chemical composition, age or distance of two stars 
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with accurate masses and radii, allowing a much more dis- 
criminating test of the physical ingredients of theoretical 
models, for example the amount of convective core over- 
shooting and mass loss, and the accuracy of the opacities 
used. Alternatively, the study of a dEB in an open cluster 
can be used to find the chemical composition of the cluster, 
its age, or its distance, avoiding the use of the technique of 
fitting theoretical isochrones to the morphology of stars in 
colour-magnitude diagrams. 

Southworth, Maxted & Smalley (2004a, hereafter Pa- 
per I) studied the early-type well-detached eclipsing binaries 
V615Per and V618Per, members of the important young 
open cluster hPersei. As the two dEBs are members of the 
same stellar group, this allowed the masses and radii of four 
stars of the same age, chemical composition and distance 
to be compared to the predictions of theoretical models, a 
situation not previously achieved. From this procedure the 
metal abundance of hPersei was found to be Z ~ 0.01, in 
conflict with previous assumptions that the cluster has an 
approximately solar chemical composition. 

Southworth, Maxted & Smalley (2004b, hereafter Pa- 
per II) studied the high-mass dEB V453 Cyg, a member of 
the sparse open cluster NGC6871. Fhe masses and radii of 
this somewhat evolved system (\ogg = 3.731,4.005) were 
determined to accuracies of 1.5%. Although the system con- 
sists of two quite dissimilar stars (with radii of 8.55 and 
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5.49 R©), current theoretical models of the Granada (Claret 
1995), Padova (Bressan et al. 1993), Geneva (Schaller et 
al. 1992) and Cambridge (Pols et al. 1998) groups showed 
an impressive agreement with the stellar properties of 
V453 Cyg. The Granada models also agreed well with the 
central condensation of the primary star derived from anal- 
ysis of the apsidal motion of the dEB. This good agreement 
suggests that we must study dEBs which are more evolved, 
of higher mass, or more dissimilar in order to assess the suc- 
cess of the different physical ingredients of stellar models. 

V621 Per is a member of the young open cluster x Persei 
and is more evolved, and composed of more dissimilar stars, 
than V453 Cyg. \ Persei is also regarded as a physical re- 
lation of h Persei (Paper I) so a full analysis of this system 
would allow the simultaneous comparison of the observed 
masses and radii of six stars (the components of V615 Per 
and V618Per, Paper I, and V621 Per, this work) with theo- 
retical predictions. 

1.1 V621 Persei 

V621 Per (Table^ was discovered to be a dEB by Krzesihski 
& Pigulski (1997, hereafter KP97) from approximately 1200 
images, through the broad-band B and V niters, of the nu- 
cleus of x Persei- The eclipses are total, last for approxi- 
mately 1.3 days, and are about 0.12 mag deep in both B and 
V . The ascending and descending branches of one eclipse 
were observed in BV on two successive nights but the only 
other observations during eclipse were 102.1 days earlier, 
and during totality, so the period could not be determined. 

The B2 giant component of V621 Per is one of the 
brightest members of \ Persei and has been studied sev- 
eral times using high-resolution optical spectroscopy to de- 
termine accurate chemical abundances. Lennon, Brown & 
Dufton (1988) derived a normal helium abundance but state 
that different lines gave different results, which they claim 
could be due to the high surface gravity used [log g = 3.6 
(cms -1 )]. Dufton et al. (1990) found a normal abundance 
of helium and various metals, but a deficiency of 0.4 dex in 
nitrogen and aluminium. These authors may also have been 
the first to note that V621 Per is a spectroscopic binary. 

Vrancken et al. (2000) derived a precise effective tem- 
perature and surface gravity of T c a = 22 500 ± 500 K and 
logg = 3.40 ± 0.05, based on the silicon ionisation balance 
and direct fitting of the H/3 and H7 absorption lines using 
non-LTE model atmosphere calculations. They also derived 
a high microturbulent velocity of 9kms _1 (or 13kms _1 
from the 11 lines) consistent with evolution away from the 
main sequence. It is notable that abundance analyses gen- 
erally find larger microturbulence velocities than those usu- 
ally assumed (see e.g., Dufton, Durrant & Durrant 1981). 
Vrancken et al. derived abundances of C, N, O, Mg, Al 
and Si comparable to other bright B stars in \ Persei, but 
abundances of the overall sample seem to be lower than 
the Sun by 0.5 ± 0.2 dex. Venn et al. (2002) derived a 
boron abundance using ultraviolet spectra taken with the 
STIS spectrograph on board the Hubble Space Telescope. 
They found a lower microturbulent velocity of 4kms _1 , as 
usual in the ultraviolet wavelength region, but a macrotur- 
bulent velocity of 20kms~ 1 . They also report a value of 
[^] = —0.16 ± 0.17dex from abundance analyses of the 
spectral lines of light metals. 



Table 1. Identifications and photometric indices for V621 Per 
from various studies. All photometric parameters refer to the com- 
bined system light (although the secondary star is much fainter 
than the primary). Most photometric quantities have been deter- 
mined many times and the quoted values have been selected as 
the most representative of all determinations. 
* Calculated from the system magnitude in the V filter, the 
adopted cluster distance modulus and reddening (see section fl. 21 
and the canonical reddening law Ay = 3.1Eg_y. 
References: (1) Argelander (1903); (2) Oosterhoff (1937); (3) 
Keller et al. (2001); (4) Slesnick et al. (2002); (5) Capilla & 
Fabregat (2002); (6) Two Micron All Sky Survey; (7) Crawford, 
Glaspey & Perry (1970); (8) Uribe et al. (2002) based on proper 
motion and position. 
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1.2 x Persei 

The open clusters \ Persei (NGC 884) and h Persei 
(NGC 869) together form the Perseus Double Cluster. The 
co-evolutionary nature of h and \ Persei has been studied 
many times since the seminal work of Oosterhoff (1937). 
The results of recent photometric studies (Marco & Bern- 
abeu 2001; Keller et al. 2001; Slesnick, Hillenbrand & Massey 
2002, Capilla & Fabregat 2002) seem to be converging to 
identical values of distance modulus (11.70 ± 0.05 mag) and 
age (logr = 7.10 ± 0.01 years), which implies that h and 
X Persei are physically related although there do appear to 
be some differences in stellar content, for example the large 
number of Be stars in \ Persei. The reddening of x Persei is 
Eb-v = 0.56 ± 0.05, but h Persei displays differential red- 
dening. These issues were discussed in detail in Paper I. 



2 OBSERVATIONS 

Spectroscopic observations were carried out in 2002 October 
using the 2.5 m Isaac Newton Telescope (INT) on La Palma. 
The 500 mm camera of the Intermediate Dispersion Spectro- 
graph (IDS) was used with a holographic 2400 Imm" 1 grat- 
ing and an EEV 4k x 2k CCD. The resulting spectra have a 



© 0000 RAS, MNRAS 000, 000-000 



Eclipsing binaries in open clusters. III. V621 Per in \ Persei 3 



signal to noise ratio of approximately 60 per pixel, measured 
in regions with no detectable spectral lines. From measure- 
ments of the full width half maximum (FWHM) of arc lines 
taken for wavelength calibration we estimate that the res- 
olution is approximately 0.2 A. The main spectral window 
chosen for observation was 4230-4500 A. This contains the 
Mgn 4481 A line which is known to be one of the best lines 
for radial velocity work for early-type stars (Andersen 1975; 
Kilian, Montenbruck & Nissen 1991). The Hei 4471 A and 
H7 (4340 A) lines are useful for determination of effective 
temperatures and spectral types for such stars. One spec- 
trum was observed at H/3 (4861 A) to provide an additional 
indicator of effective temperature. 

Data reduction was undertaken using optimal extrac- 
tion as implemented in the software tools PAMELA and 
molly 1 (Marsh 1989). 



3 SPECTROSCOPIC ORBIT 

The INT spectra contain identifiable spectral lines from only 
the primary star. Radial velocities were derived from the 
spectra by cross-correlation with a synthetic template spec- 
trum, using the XCOR routine in molly. Several template 
spectra were investigated and the resulting radial velocities 
were found to be insensitive to the choice of template. 

As our spectroscopic observations cover less than the 
full orbital period of V621 Per, and the orbit is eccentric, 
they cannot provide a unique value of the period. To solve 
this problem, literature radial velocities were taken from Liu, 
Janes & Bania (1989, 1991) and Venn et al. (2002). Addi- 
tional high-resolution spectra were generously made avail- 
able by Dr. P. Dufton and Dr. D. Lennon. These were orig- 
inally observed to determine chemical abundances (Dufton 
et al. 1990; Vrancken et al. 2000) so the wavelength cali- 
brations may only be accurate to a few kms -1 (P. Dufton, 
private communication). Radial velocities were derived by 
fitting Gaussian functions to strong spectral lines, predom- 
inantly from the He 1 and 11 ions, and excellent agreement 
was found between different lines in the same spectra. 

Using the photometric constraint that the orbital period 
of V621Per must be a submultiple of 102.1 days (KP97), 
the possible periods were investigated by analysing the ra- 
dial velocities with SBOP 2 and the KP97 light curves with 
ebop (SectionJSJ. Only a period around 25.5 days can pro- 
vide a good fit to all the data. The photometric and spec- 
troscopic data were fitted simultaneously by requiring the 
spectroscopically-derived orbital period to correctly predict 
the phase at which the primary eclipse occurs. The resulting 
orbital ephemeris is 

T pcri = HJD 2 452 565.150(97) +25.53018(20) x E 

where T per i is the time of periastron passage and the num- 
bers in parentheses refer to the uncertainty in the last digit 
of the quantity. The midpoints of the primary and secondary 
eclipses occur at phases 0.67 and 0.06, respectively. 

The final spectroscopic orbit was calculated using the 

1 PAMELA and MOLLY were written by Dr. Tom Marsh and are 
found at http://www.astro.soton.ac.uk/~trm/software.html 

2 Spectroscopic Binary Orbit Program written by Dr. Paul B. 
Etzel (http : //mintaka. sdsu. edu/f aculty/etzel/). 



Table 2. Radial velocity observations of V621 Per and the O — C 
values with respect to the final spectroscopic orbit. 
References: (1) Liu, Janes &: Bania (1989); (2) Liu, Janes & 
Bania (1991); (3) measured from spectra generously provided by 
Dr. P. Dufton; (4) Venn et al. (2002); (5) This work. 
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radial velocities derived from the INT spectra and fixing the 
orbital period at the value given above. The orbit is plotted 
in Fig.0 the radial velocities and observed minus calculated 
(O — C) values are given in Table [5] and the parameters of 
the orbit are given in Table |3] 

The projected rotational velocity of the primary com- 
ponent of V621 Per was found by fitting Gaussian functions 
to the Sim 4575 A spectral line singlet. Using the orbital 
inclination found in Section |S] the equatorial rotational ve- 
locity is Voq = 32.2 ± 1.2 kms - where the quoted error is 
the 1 a error of the individual values. 

The INT spectra are single- lined in character, and in 
experiments with the two-dimensional cross-correlation al- 
gorithm TODCOR (Zucker & Mazeh 1994) and with spec- 
tral disentangling (Simon & Sturm 1994) we were unable 
to detect any signal from the secondary star. By simulat- 
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Figure 1. Spectroscopic orbit for V621 Per. Filled circles denote 
radial velocities derived from the INT spectra and open circles 
show radial velocities obtained from other sources. The systemic 
velocity is indicated by a dotted line. 



Table 3. Parameters of the spectroscopic orbit derived for 
V621Per. 
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ing the spectrum of the secondary star with a rotationally 
broadened primary spectrum we have constructed several 
trial composite spectra. From analysis of these using cross- 
correlation, we estimate that we would have detected sec- 
ondary spectral lines if it contributed more than 5% of the 
total light, for a rotational velocity of 50 kms -1 . If it rotates 
faster than this, or has a spectrum very different to that of 
the primary star, then the detection threshold will increase. 



4 DETERMINATION OF EFFECTIVE 

TEMPERATURE AND SURFACE GRAVITY 

Lennon, Brown & Dufton (1988) found the atmospheric pa- 
rameters of V621 Per to be T cff = 21 500 K and logg = 3.6, 
from several different Stromgren photometric calibrations 
and fitting Balmer lines with synthetic profiles. Dufton et 
al. (1990) found T eS = 21 700 K and log 5 = 3.6 using a 
similar method, but found T e g = 23 000 K from the silicon 
ionization equilibrium, for which the corresponding \ogg is 
3.7. 

Vrancken et al. (2000) derived T cff = 22 500 ± 500 K 
and log g = 3.40 ± 0.05 from the silicon ionization equilib- 
rium and fitting Balmer lines with synthetic profiles. These 
atmospheric parameters were adopted by Venn et al. (2002) . 

Using Stromgren uvbyf3 data taken from Crawford, 
Glaspey & Perry (1970), and the calibration of Moon & 



Dworetsky (1985), we derive T cS = 21700 ± 800 K and 
logg = 3.69 ± 0.07 (where the uncertainty is a formal error 
of the fit). The Stromgren photometry of Marco & Bernabeu 
(2001) gives T cff = 19 300 ± 800 K and logg = 3.56 ± 0.07, 
and of Capilla & Fabregat gives T cB = 20 900 ± 800 K and 
logg = 3.36 ± 0.07. The Crawford et al. (1970) photometry 
should be preferred as the filters are closest to the filters used 
to define the Stromgren uvby and Crawford /3 systems, and 
because photoelectric uvby/3 photometry has been shown to 
be superior to CCD uvby/3 photometry (see e.g., Mermilliod 
& Paunzen 2003). 

Using Geneva photometry from Rufener (1976) and the 
calibration of Kunzli et al. (1997) we find T cff = 22 230 ± 
250 K and a high log g value of 3.97±0.18. Use of the Geneva 
photometry of Waelkens et al. (1990) gives T cS = 23 200 ± 
420 K and a low log 5 value of 3.31 ± 0.26. 

We have assumed a surface gravity value of log g = 3.6, 
which agrees well with Dufton et al. (1990), the calibra- 
tion results using the Crawford et al. (1970) data, and with 
our spectroscopic and photometric analyses (see Section |SJ. 
We have fitted our H7 and H/3 spectra with synthetic pro- 
files calculated using uclsyn (Smith 1992, Smalley et al. 
2001; see Paper I for further references) and Kurucz (1993) 
ATLAS9 model atmospheres. The spectra were rotationally 
broadened as necessary and instrumental broadening was 
applied with FWHM = 0.2 A to match the resolution of the 
observations. For log g = 3.60 we find T cff = 22 500 ± 500 K, 
in agreement with Vrancken et al. (2000). 



5 LIGHT CURVE ANALYSIS 

We have analysed the BV light curves of KP97 using the 
simple and efficient NDE eclipsing binary model (Nelson 

6 Davis 1972) as implemented in the light curve analysis 
program ebop 3 (Popper & Etzel 1981). In this model the 
stellar shapes are approximated by biaxial ellipsoids. The 
orbital eccentricity and longitude of periastron were fixed 
at the spectroscopic values, initial filter-specific linear limb 
darkening coefficients of 0.30 (primary star) and 0.25 (sec- 
ondary star) were taken from van Hamme (1993) and the 
gravity darkening exponents /3i were fixed at 1.0 (Claret 
1998). Changes in the limb darkening and gravity brighten- 
ing values for the secondary star have a negligible effect on 
the photometric solutions because this star contributes very 
little of the light of the system. Contaminating 'third' light 
was fixed at zero, as solutions in which it was a free param- 
eter were not significantly different from solutions with no 
third light. 

Initial light curve solutions converged to an orbital in- 
clination, i, of 90°, but values of i from about 88 to 90° 
fit the observations almost equally well. Solutions in which 
the surface brightness ratio, J, was freely adjusted towards 
the best fit generally converged to a value of J below zero, 
which is unphysical. We therefore present separate solutions 
(Tabled for the B and V light curves in which orbital incli- 
nation is fixed at i = 88, 89 and 90° and the surface bright- 
ness ratio is fixed at J = 0.0, 0.5 and 1.0. The light curves 

3 Eclipsing Binary Orbit Program written by Dr. Paul B. Etzel 
(http://mintaka.sdsu.edu/faculty/etzel/). 
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Table 4. Results of the light curve analysis of V621 Per for several different (fixed) values of the surface brightness ratio and orbital 
inclination. The final entry gives the adopted values and uncertainties of the parameters (see text for discussion). The uncertainties are 
confidence intervals, not 1 a errors. 
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Figure 2. The KP97 B and V light curves of V621 Per around the primary eclipse, phased using the spectroscopic ephemeris, with the 
best-fitting EBOP model light curves. The B light curve, offset by —0.1 mag for clarity, is shown using open circles and the V light curve 
is shown using filled circles. The best-fitting curves were generated with J = 0.25 and i = 89.0°. 
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Figure 3. Results of the Monte Carlo analysis for J = 0.25 
and i = 89.0°. The limb darkening coefficients, «a and ub, were 
chosen randomly on a flat distribution between u — 0.05 and u + 
0.05 for each synthetic light curve, and fixed during solution of 
the light curve. 



are of insufficient quality to solve for the limb darkening 
coefficients so these have been fixed during solution. 

Robust errors were estimated using Monto Carlo sim- 
ulations (see Paper II for a detailed description of the 
method). In this procedure, a best fitting model is found 
for the observed light curve. Ten thousand synthetic light 
curves are then generated by evaluating the model bright- 
ness at each phase of observation, and adding observational 
noise of the same magnitude as that present in the observed 
light curve. Each synthetic light curve is fitted, and the re- 
sults of the fits allow a robust estimation of the uncertainties 
and of the parameter correlations at the point of best fit. 
This method was modified to explicitly include uncertain- 
ties due to the use of assumed limb darkening coefficients by 
fixing them at random values on a flat distribution within 
±0.05 of the original value. 

As no trace of the secondary star was found in the 
observed spectra, the B-filter light ratio must be 0.05 or 
less. The maximum light ratio in the V filter will be slightly 
greater than this as the secondary star is expected to have a 



lower effective temperature than the primary star. For sim- 
plicity, we have adopted a maximum light ratio of 0.05 for 
both the B and V light curves. We have therefore calcu- 
lated best-estimate parameters by evaluating the ranges of 
possible parameter values in the two light curves and then 
averaging the midpoints of the ranges for the two light curves 
(Table [IJ. The quoted uncertainties are confidence intervals 
which encompass the range of possible values for each pa- 
rameter, so are not 1 a errors. This procedure is simple but 
is quite adequate considering the nature of the observations 
analysed here. 

Fig. shows the observed light curves and the best- 
fitting models with J = 0.25 and i = 89°. Fig. [3] represents 
the relation between different parameters of the fit to the 
V light curve. As V621 Per exhibits total eclipses, the radii 
of the two stars are only weakly correlated. Changes in the 
limb darkening coefficients used do affect the derived radii of 
both stars, but this effect is quite small and easily quantified. 



6 ABSOLUTE DIMENSIONS AND 

COMPARISON WITH STELLAR MODELS 

Although the absolute masses and radii of the component 
stars of V62f Per cannot be found directly, the mass function 
and fractional radii (the stellar radii expressed as a fraction 
of the semi-major axis of the orbit) are accurately known. 
This allows us to empirically determine the surface gravity 
of the secondary star despite not knowing its actual mass or 
radius. 

From consideration of the definitions of the mass func- 
tion, fractional radius and surface gravity, and from Kepler's 
third law, it can be shown that the surface gravity of the 
components of a single-lined dEB can be given by 



9-n = 



2ttN 4 / 3 [G/(M)] 1/3 
P ) r„ sin i 



(M n \ 
\M B ) 



(1) 



where P is the orbital period, f(M) is the mass function, 
r n and M n are the fractional radius and absolute mass of 
component n (n = A, B) and Mb is the mass of the spec- 
troscopically unseen star. In the usual astrophysical units of 
solar masses, days and gravity in cms -1 , the surface gravi- 
ties of the two stars are therefore given by 



loggB = 3.18987 + 



log/(M) 41ogP 



log(r B sini) (2) 



logfifA = 3.18987+ 



log/(Af) 



41ogP . 2 . .. 

log(r A smi)-logq(3) 



3 3 

where q — is the mass ratio. Equation [5] contains only 
known quantities so, despite not knowing the mass or radius 
of the secondary component of V621 Per we can empirically 
calculate its surface gravity to be log gB = 4.244 ± 0.054. 

Alternatively, it is possible to use V621 Per's member- 
ship of the open cluster xPersei to infer the properties of 
the primary star. The absolute magnitudes of the V621 Per 
system, found from the apparent magnitudes, the distance 
modulus and reddening of the cluster (Table Q and Sec- 
tion ll.2l and the reddening laws Av = 3AEb-v and Ak = 
0.38E B -v (Moro & Munari 2000), are M v = -4.04 ± 0.16 
and Mk = —3.20 + 0.06. Adopting bolometric corrections of 
-2.20 + 0.05 and -2.93 + 0.08 (Bessell, Castelli & Plez 1998) 
gives absolute bolometric magnitudes of —6.24 ± 0.17 and 
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Figure 5. Hertzsprung-Russcll diagram showing the luminos- 
ity and effective temperature derived for V621 Per. The Granada 
evolutionary model preditions are plotted for for masses of 10.0, 
12.8 and 15.8 Mq, and for metal abundances of Z = 0.01 (dashed 
lines) and 0.02 (dotted lines). The zero age main sequences are 
plotted using the same line styles and filled circles denote predic- 
tions for ages of 12.3, 12.6 and 12.9 Myr. 

—6.13 ± 0.10 for the V and K-RlteT data respectively. The 
two values are in good agreement but the 7v-filter value is 
more accurate because it is less affected by the uncertainty 
in Eb-v- We will adopt the A'-filter value as it is more ac- 
curate, and the 2MASS apparent magnitudes are known to 
be very reliable. Adopting a solar absolute bolometric mag- 
nitude of 4.74 (Bessell et al. 1998), which is consistent with 
the adopted bolometric corrections, gives a luminosity of 
log ^ = 4.348±0.039. This gives a radius of 9.9±0.7 R for 
the primary star. If we assume that the secondary star's con- 
tribution to the total light of the system is 5%, this will cause 
the primary radius to be overestimated by about 0.25 Rq, 
which is negligible at this level of accuracy. 

6.1 Comparison with stellar models 

The mass ratio of V621 Per can in principle be found from 
equation[3]as log^A has been determined by several authors 
(see Sectional . However, the mass ratio has a very sensitive 
dependence on the value of log <?a as it is effectively calcu- 
lated from the difference of two numbers of similar size. An 



alternative approach is to evaluate the mass ratio, and hence 
the absolute masses and radii of the two stars, for several 
different values of log gA and compare the possibilities with 
the predictions of stellar models. Substituting the mass ratio 
into the definition of the mass function, we can derive the 
absolute masses of the two stars using 

M A =^^ (4) 

Mb = qM A (5) 

The absolute radii of the stars can then be found from the 
masses and the surface gravities. 

We have used the equations above to determine the 
absolute masses and radii of both components of V621 Per 
for several assumed values of log<?A between 3.40 and 3.70. 
These have been compared to the predictions of the Granada 
stellar evolutionary models (Claret 1995; Claret & Gimenez 
1995, 1998) for ages close to the age of the x Persei open 
cluster, 12.6 ± 0.3 Myr (Paper I). 

Fig.^Jshows possible values of the absolute masses and 
radii of the components of V621 Per compared to predic- 
tions of the Granada evolutionary models for ages around 
12.6 Myr. Also shown (by a shaded area) is the range of pos- 
sible values of the primary radius, derived from the known 
distance and apparent magnitude of the dEB. Whilst all 
three possible metal abundances can fit the two stars for 
loggA ~ 3.55, the Z — 0.01 predictions provide the best 
fit to the predicted properties of the secondary component. 
This diagram suggests that logpA is probably between 3.55 
and 3.60. If the best fit is sought for metal abundances of 
Z = 0.02 and 0.004, ages of roughly 5 and 40 Myr, respec- 
tively, are found. 

The primary component of V621 Per has been plot- 
ted in the Hertzsprung-Russell diagram (Fig. and com- 
pared with Granada theoretical model predictions for masses 
of 10.0, 12.8 and 15.8 Mq and for metal abundances of 
Z — 0.01 and 0.02. The \ Persei open cluster has been found 
to have an age of 12.6 ± 0.3 Myr (Paper I) and this age has 
been indicated on the evolutionary track for for each model 
mass. The position of the primary component of V621 Per 
in these diagrams suggests that its mass is a little below 
12.8 Mq, also consistent with the form of the mass-radius di- 
agram (Fig.^J. However, its age derived by comparison with 
evolutionary models is somewhat greater than the 12.6 Myr 
expected due to its membership of \ Persei, and the discrep- 
ancy is larger compared to the Z = 0.01 model predictions 
than for the Z = 0.02 predictions. 

The age of 12.6 ± 0.3 Myr for \ Persei was derived, from 
comparison between photometric observations of the cluster 
and the predictions of theoretical models, by researchers who 
assumed that Z — 0.02 (Paper I). As the metal abundance 
of x Persei has been found to be Z = 0.01 (Paper I), this age 
may have a systematic error. Therefore the age discrepancy 
found here is of only minor significance, but deserves inves- 
tigating when more accurate parameters are determined for 
V621 Per. Also, the amount of overshooting present in stel- 
lar models is known to significantly change the predicted 
ages of giant stars (e.g., Schroder & Eggleton 1996) and the 
inclusion of rotation also affects the main sequence lifetime 
of high- mass stars (Maeder & Meynet 2000). 
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Table 5. Absolute masses and radii of the components of V621 Per calculated using different values of logg^- The luminosity is that of 
the primary component only as we do not know the effective temperature of the secondary star. 



loggA 

(cms ) 


Mass ratio 


Primary mass 
(M©) 


Secondary mass 
(M ) 


Primary radius 

(Re) 


Secondary radius 

(R©) 


Luminosity 
(logL/L ) 


3.40 


0.676 ± 0.049 


5.6 ± 1.0 


3.8 ± 0.4 


7.8 ± 0.7 


2.4 dz 0.2 


4.152 ± 0.057 


3.45 


0.602 ± 0.044 


7.3 ± 1.3 


4.4 ± 0.4 


8.4 ± 0.7 


2.6 ± 0.2 


4.213 ± 0.058 


3.50 


0.537 ± 0.039 


9.4 ± 1.8 


5.1 ± 0.5 


9.0 ± 0.8 


2.8 dz 0.2 


4.277 dz 0.059 


3.55 


0.478 ± 0.035 


12.3 ± 2.4 


5.9 dz 0.6 


9.8 ± 0.9 


3.0 dz 0.2 


4.343 dz 0.061 


3.60 


0.426 ± 0.031 


16.2 dz 3.2 


6.9 ± 0.8 


10.6 ± 1.0 


3.3 ± 0.3 


4.412 ± 0.062 


3.65 


0.380 ± 0.028 


21.4 dz 4.3 


8.1 ± 0.9 


11.5 ± 1.1 


3.6 ± 0.3 


4.483 dz 0.063 


3.70 


0.339 ± 0.025 


28.5 ± 5.9 


9.6 dz 1.2 


12.5 ± 1.2 


3.9 dz 0.3 


4.557 dz 0.064 




0.6 0.8 1.0 1.2 1.4 

log(M/M ) 



Figure 4. The logarithmic mass— radius plot for the two components of V621 Per. The possible combinations of mass and radius for 
each star, for different values of loggAi are plotted with errorbars connected by solid lines. Numbers on the diagram indicate the value of 
loggA used to calculate the adjacent datapoint. The Granada stellar model preditions are plotted for an age of 11.0 Myr for Z = 0.004 
(dash-dotted line), Z = 0.01 (dashed line) and Z = 0.02 (dotted line). The 11.5 Myr Z = 0.01 predictions are plotted using a dashed line 
with the symbols representing the points of model evaluation omitted. Curves have been calculated using a cubic spline interpolation. 
The radius of the primary star, calculated from its known distance and apparent magnitude, has been shown using light shading to 
indicate the range of possible values. 



From comparison with the Granada evolutionary mod- 
els, the surface gravity of the primary component of 
V621 Per is approximately 3.55. This conclusion is valid for 
metal abundances of Z = 0.01 and 0.02. From Tabled the 
masses and radii of V621 Per corresponding to log g\ = 3.55 
are about 12 and 6Mq, and 10 and 3Rq, for primary and 
secondary star respectively. This means that the primary 



star is near the age at which it passes through the 'blue 
loop' evolutionary stage (the point at which core hydro- 
gen exhaustion causes the effective temperature and surface 
gravity to rise temporarily), so accurate masses and radii 
for it would provide extremely good tests of the predictions 
of theoretical models. The Granada stellar models predict 
a luminosity ratio of about 0.05 for the inferred properties 
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of V621 Per. The light ratio in the blue will be smaller than 
the overall luminosity ratio because the secondary star has 
a lower effective temperature than the primary star. This 
suggests that the quality of our spectroscopic observations 
was almost sufficient to detect the secondary star. Accurate 
velocities for both stars should be measurable on spectra 
of a high signal to noise ratio, depending on the rotational 
velocity of the secondary star. 

6.2 Membership of the open cluster ^Persei 

V621 Per is situated, on the sky, in the centre of the \ Persei 
open cluster. It also appears in the correct place on the 
colour-magnitude diagram of the cluster (see Section II. 21 
and has the correct proper-motion (Uribe et al. 2002) for 
cluster membership. The systemic velocity of the dEB, 
—44.5 ± 0.4 km s^ 1 , is consistent with the measured cluster 
systemic velocities of Oosterhoff (1937), Bidelman (1943), 
Hron (1987), Liu et al. (1989, 1991) and Chen, Hou & Wang 
(2003). In Paper I we measured the systemic velocity of the 
hPersei cluster to be —44.2 ± 0.3 kms -1 , indicating that h 
and x Persei have the same systemic velocities, which is con- 
sistent with them having a common origin. 



7 DISCUSSION 

V621 Persei is a detached eclipsing binary in the young open 
cluster x Persei, composed of a bright B2 giant star and 
an unseen main sequence secondary star. From blue-band 
spectroscopic data and radial velocities taken from the liter- 
ature, we have derived an orbital period of 25.5302 days and 
a mass function of f(M) = 0.617 ± 0.012 M Q . The discovery 
light curve of KP97 shows that the system exhibits a total 
primary eclipse lasting around 1.3 days and about 0.12 mag 
deep in B and V . No data exist around phase 0.06, where 
the secondary eclipse is expected to occur. The secondary 
eclipse may be up to about 0.06 mag deep. The light curves 
have been solved using ebop and Monte Carlo simulations 
to find robust uncertainties, and accurate fractional radii 
have been determined. The surface gravity of the secondary 
component has been found to be logpe = 4.244 ± 0.054. 

Using the data above, possible values of the absolute 
masses and radii of the two stars were calculated by as- 
suming different values of the primary surface gravity. A 
comparison in the mass-radius diagram of these possible 
values with theoretical predictions from the Granada stel- 
lar evolutionary models suggests that log^A ~ 3.55. This 
surface gravity value agrees well with the values determined 
by Lennon et al. (1988) and Dufton et al. (1990) by fitting 
observed Balmer line profiles with synthetic spectra. The 
luminosity of V621 Per has been derived from the known 
distance of the x Persei cluster and the apparent magnitude 
of the dEB. This has been used to place the primary star 
in the Hertzsprung-Russell diagram, and a comparison with 
the Granada evolutionary models confirms that its mass is 
roughly 12 Mq, although a small discrepancy exists between 
the inferred age of V621 Per and the age of x Persei. 

The value of log g\ leads to masses of approximately 12 
and 6 Mq and radii of 10 and 3 R© for the components of 
the dEB. This conclusion is not strongly dependent on use 
of the Granada stellar models; predictions of the Geneva, 



Padova and Cambridge models are close to those of the 
Granada models (Paper II). V621 Persei is a potentially im- 
portant object for the information it holds about the evo- 
lution of high-mass stars. The expected luminosity ratio of 
the system, about 0.05, suggests that spectral lines of the 
secondary component should be detectable in high signal-to- 
noise spectra. Better light curves will be needed for detailed 
studies of the properties of V621 Per, but the long period 
and lengthy eclipses mean that a large amount of telescope 
time will be required. In particular, the secondary eclipse 
is expected to occur around phase 0.06 and may be up to 
0.06 mag deep. Observations of the light variation through 
secondary eclipse will be needed to provide a definitive study 
of the system. The absolute dimensions of the primary star, 
a B2 giant which is close to the blue loop stage of stellar 
evolution, could provide a good test of the success of stellar 
evolutionary models and of the amount of convective core 
overshooting which occurs in stars. 
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